Quantum mechanical potentials related to the prime numbers and Riemann zeros 
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Prime numbers are the building blocks of our arithmetic, however, their distribution still poses 
fundamental questions. Bernhard Riemann showed that the distribution of primes could be given 
explicitly if one knew the distribution of the non-trivial zeros of the Riemann ^(s) function. Ac- 
cording to the Hilbert-Polya conjecture there exists a Hermitean operator of which the eigenvalues 
coincide with the real part of the non-trivial zeros of C,{s). This idea encourages physicists to exam- 
ine the properties of such possible operators, and they have found interesting connections between 
the distribution of zeros and the distribution of energy eigenvalues of quantum systems. We apply 
the Marchenko approach to construct potentials with energy eigenvalues equal to the prime numbers 
and to the zeros of the C,{s) function. We demonstrate the multifractal nature of these potentials 
by measuring the Renyi dimension of their graphs. Our results offer hope for further analytical 
progress. 

PACS numbers: 02.10.De, 02.30.Zz, 03.65.Ge 
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I. PRIMES, ZEROS AND QUANTA 

The prime numbers are the building blocks for the pos- 
itive integers, since the fundamental theorem of arith- 
metic states that every positive integer can be written as 
a product of primes, and this product is unique up to a 
rearrangement of the factors. Additionally, not only does 
the product of prime numbers have this remarkable prop- 
erty, but also their sum. Only a decade ago O. Ramare 
proved [ij that any positive integer can be written as a 
sum of no more than six prime numbers, but the Gold- 
bach conjecture [31, whether every number is expressible 
as a sum of two prime numbers, remains unproven. 

Based on empirical evidence, many mathematicians 
conjectured that the prime counting function, t:{x) = 
I {p I p is prime and p < asymptotically behaves as 
the logarithmic integral Li(a;). Hadamard 0] and de la 
Vallee Poussin [4] independently gave the rigorous proof 
for this statement. Riemann derived the following exact 
formula @ 



it{x) = lim R{x) 




(1) 



where R{x) is the so-called Riemann function defined as 
= X]m=i /i(m) denotes the 
Mobius function. The sum in ([1]) is extended over all non- 
trivial zeros, p, of the Riemann C(s) function, counted 
with their multiplicities. The latter function ^(s) is de- 
fined by the infinite series, = ^O'' s > 1, 
and, otherwise, by its analytic continuation over the com- 
plex s plane. 
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Exploring the locations of the zeros of C,{^)i Riemann 
made his famous conjecture: all the non-trivial zeros lie 
on the s = i + it is real) critical line. Proving or dis- 
proving the Riemann-hypothesis remains the most tan- 
talising challenge in number theory since David Hilbert 
nominated it in 1900 Q as the eighth problem on his 
famous list of compelling problems in mathematics [27] . 

The connection between the Riemann-hypothesis and 
physics seems to date back to the early years of quantum 
mechanics. According to the Hilbert-Polya conjecture, 
the zeros of C,{s) can be the spectrum of an operator, 
O = -iX + iH, where H is self-adjoint. This operator Ti, 
might have a physical interpretation as a Hamiltonian of 
a physical system and, therefore, the key to the proof of 
the Riemann Hypothesis may have been coded in physics. 

The analogy between the properties of <^(s) and the en- 
ergy eigenvalues of a quantum mechanical system provide 
us some information about the form of a possible opera- 
tor Ti Q. One of these similarities, the comparison of the 
number of ^(s) zeros and the number of energy eigenval- 
ues below a threshold, suggests that the physical system 
is quasi-one-dimensional. This link is further strength- 
ened by checking different statistics of the zeros, such 
as the nearest-neighbour spacings, the n-correlations be- 
tween the zeros, etc. Montgomery showed that these dis- 
tributions are all in good agreement with the Gaussian 
unitary ensemble of random matrix theory ,8] . 

In this work we utilize an inverse scattering formalism 
and construct potentials of which the energy eigenval- 
ues are the zeros of the Riemann C(s) function. We also 
consider the problem when the eigenvalues are taken to 
be the prime numbers themselves. In Sec. |TT] we intro- 
duce our numerical method and give evidence that it is 
capable of generating potentials from sets of discrete en- 
ergy eigenvalues, such as a finite set of C, zeros or a finite 
set of prime numbers. We calculate the Renyi dimension 
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[Pi for these potentials. The results anticipate that these 
potentials have multi-fractal nature. In section Hill we at- 
tempt to further clarify why previous studies presented 
contradictory results for the fractal dimension. 



II. INVERSE SCATTERING FORMALISM 

Provided the Hilbert-Polya conjecture is true, the nat- 
ural and plausible approach to finding operator Ti would 
be to approximate it from a finite number of eigenvalues. 
We will follow this path assuming the existence of a local 
potential V whose spectrum is related to the zeros of the 
Riemann ^(s) function or, later, to the prime numbers. 

The Marchenko approach aims to reconstruct a sym- 
metric potential using the characteristics of both the 
bound states (energy eigenvalues and normalisation con- 
stants) and the scattering states (reflection coefficient 
at all energies). The question of the existence and 
uniqueness of any solution obtained from the inversion 
procedures is delicate, although if one assumes a one- 
dimensional, symmetric potential the complete set of 
eigenvalues uniquely determines the potential [lo| . Dif- 
ferent, but mathematically equivalent, methods exist 
for reconstructing the scattering potential in a one- 
dimensional quantum mechanical problem. 

For a given set of energy eigenvalues and reflection co- 
efficient the quantum potential can be calculated from 



V{x)^-2—K{x,x% 
ax 



(2) 



where the order of operations is important. One should 
localise x' first and then perform the differentiation. The 
function K{x,x') is the solution of the Marchenko inte- 
gral equation {x' > x) 

1*00 

K{x,x')+Ko{x+x')+ K{x,s)Ko{s + x')ds = (3) 



and the kernel K(){z) is determined by the spectral pa- 
rameters 



Ko{z) 



-, poo N 

2^ 



(4) 



The input data are the reflection coefficient, R{k), at 
energy E — h'^k^/2m, and the normalisation constant, 
c„, of the nth bound state related to the discrete energy 
eigenvalue E„ = — ?i^«;^/2m. 

In general, scattering states contribute to the kernel, 
and therefore a family of potentials can be associated 
with a given set of energy eigenvalues. However, accord- 
ing to our physical picture there are no scattering states 
in our case. The potential is theoretically inflnitely deep, 
because the set of prime numbers and the set of ({s) ze- 
ros are infinite sets with no upper bound. Henceforward 
we take R{k) = 0. 

In the case of the reflectionless potential, the scattering 
states in equation ^ do not contribute to the kernel 



and, therefore, lead to a separable and exactly solvable 
integral equation. The potential can then be obtained 
from the following formula [13] 



F(a;) ^-2-^ln(det (I-hC)) 
where I denotes identity matrix, and 



2:r = n 



TV 



m—1 



(5a) 

(5b) 
(5c) 



Choosing a set of eigenvalues, {kh}, one may calculate 
the corresponding normalisation constants, {c„} and ma- 
trix C from equations (|5bl -c). Substituting this matrix 
into ([5a|) . the desired potential can be calculated. 

We note here briefly, that using the matrix identity, 
ln(det(M)) = Tr(ln(M)), and the power series expan- 
sion of In (1 + a;), one can symbolically derive the follow- 
ing expression for the potential 



V{x) 



2V(^Tr 



(6) 



where Tr (M) denotes the trace of matrix M. 

The formulae (jSal -c) above form the basis of our cal- 
culations. Although these expressions may seem simple, 
the accurate numerical evaluation of the determinant can 
prove challenging, particularly as the number of eigenval- 
ues is increased. 

The inversion technique in its present form is numer- 
ically not convenient for more than about five hundred 
eigenvalues, for two reasons: firstly the matrices involved 
are dense, and secondly, the numerical precision required 
is demanding, since for medium values of x one has to cal- 
culate the exponential functions in (|5bp very accurately 
to have precise cancellation. The transformation of the 
formulae (j5al -c) into a numerically more tractable form is 
under investigation. 

For this reason, for large sets of eigenvalues, we used 
the dressing-transformation [T^] to calculate the poten- 
tial. We have checked numerically, up to three hundred 
energy eigenvalues, that Marchenko's inversion method 
and the dressing-transformation give identical results 
within numerical accuracy. The strength of our approach 
is in the explicit formulae for the potential construction. 
Using the dressing-transformation one has to recursively 
solve ordinary differential equations, since the potential 
is built up by incorporating the energy eigenvalues one 
by one, so in every step the solution of the previous step 
is used. Therefore, the applicability of this method to 
gain general and analytical results is limited. Contrary, 
in our method all quantities are expressed in terms of 
the input parameters, viz. the set of energy eigenvalues, 
offering some hope of analytic progress. 
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A. Reconstruction of well-known potentials 

To illustrate the method, we reconstruct well-known 
potentials from their spectra, using the triangular and 
harmonic oscillator potentials. Later the extended Nu- 
merov method [l^ is used to calculate the energy eigen- 
values of the inversion potentials. This also serves to 
check the validity of the potentials obtained using the 
Marchenko approach. 



TABLE I: Comparison of the known energy eigenvalues, eo,n, 
for the harmonic and triangular reference potentials with the 
energy eigenvalues, e„, of the appropriate inversion potentials 
calculated for the first one hundred eigenvalues. The table 
shows the first and the last five eigenvalues. 

Harmonic potential Triangular potential 



n eo,i 



£0,1 



f 
2 
3 
4 
5 

96 
97 
98 
99 

too 



f 
3 
5 
7 
9 

f9f 
f93 
f95 
f97 
f99 



f.00f923 
3.000020 
5.000944 
7.00003f 
9.000696 

f90.997520 
f92.99664f 
f94.99602f 
f97.0f5433 
f98.984293 



f.0f8793 f.0f5439 

2.338f07 2.338100 

3.248198 3.247152 

4.087949 4.087942 

4.820099 4.819397 



36.995074 
37.252699 
37.509795 
37.765659 
38.021009 



36.995066 
37.252615 
37.509785 
37.765580 
38.021020 



We build up the potentials using a finite number of 
eigenvalues, starting with one then two, five and finally 
one hundred eigenvalues from the bottom of the known 
spectrum. As one may expect, incorporating more and 
more eigenvalues into the method results in the inversion 
potential becoming more and more accurate and repro- 
ducing the spectrum faithfully. This tendency is clearly 
captured in Figure ([1]) for the triangular potential and 
in Figure ([2]) for the harmonic potential. Furthermore, 
the inversion potentials reach their asymptotic value ex- 
ponentially 12] and this asymptote lies between the last 
energy eigenvalue used for the inversion and the next 
eigenvalue. Although the discrepancy between the exact 
and inversion potentials becomes larger toward the edge 
of the inversion potential, the energy eigenvalues are still 
correctly reproduced (see Table |I| with tolerable errors. 



B. Inversion potential for prime numbers 

Using semi-classical arguments, one may show that for 
a one-dimensional potential the energy eigenvalues can- 
not increase more rapidly than quadratically, i.e. e„ ^ 

Intuitively this may be seen by noting that in the case 
of the triangular attractive potential e„ scales as r?!^ ^ 
while for the harmonic oscillator e„ goes as n, and, as a 
limiting case, for the infinite-box potential, varies as 
r? . A corollary of the Hadamard-Poussin theorem [Tsj 



TABLE II: The energy eigenvalues of the potentials derived 
for the prime numbers and the zeros of the (^(s) function. For 
both the primes and zeros of the Riemann ("(s) function, the 
first column comprises the exact eigenvalues, eo,n, the second 
one contains tsc.n for the semi-classical potential, and the last 
incorporates the energy eigenvalues e„ of the appropriate in- 
version potential. 

Prime numbers Riemann ("(s) zeros 



n eo,; 



eo,; 



2 
3 
5 
7 
11 



0.6895 
2.5316 
5.0674 
7.9717 
11.1201 



1.6387 
3.0005 
4.7052 
7.0006 
10.702 



14.1347 
21.0220 
25.0109 
30.4249 
32.9351 



13.4690 
23.2274 
29.8790 
36.0644 
41.4187 



13.0302 
21.0208 
24.7026 
30.4234 
32.8091 



96 503 513.8440 503.0008 

97 509 520.3027 508.7052 

98 521 526.7728 520.9981 

99 523 533.2544 522.9371 
100 541 539.7472 540.9843 



229.3374 284.3914 229.3354 

231.2502 287.1530 231.2259 

231.9872 289.9657 231.9865 

233.6934 292.8088 233.6322 

236.5242 295.7020 236.5215 



is that the nth prime number is approximately nln(n), 
which is clearly less than n^. We cannot, therefore, 
rule out the existence of a quantum mechanical poten- 
tial which has prime numbers as energy eigenvalues. 

We now turn to the construction of a semi-classical po- 
tential for which the first n energy eigenvalues coincide 
with the first n prime numbers. There is no theoreti- 
cal limit on the number of incorporated prime numbers, 
although numerically the calculation becomes quite cum- 
bersome. 

Using the Wentzel-Kramers-Brillouin semi-classical 
quantisation formula J^] and the leading terms of the 
prime number counting function, 7r(i?) k, R{E), one may 
derive Yf\ the following implicit equation for a poten- 
tial of which the eigenvalues are approximately the prime 
numbers {Eq > 1): 



m— 1 



TO 



E- 



\u{E)^/V-E 



■.dE 



(7) 



where E^ denotes the reference energy-level. Due to the 
density of prime numbers, p{E) w \n{E), this reference 
energy cannot be less than 1. Even though the inte- 
gral, for general m, cannot be expressed using elemen- 
tary functions, one may bound the integral from below 
and from above such that (x ^ 1) 



' In^ (x) < V{x) - Eo<x^ In^ {xln{x)). 



(8) 



In Figure [3] we plot the inversion potential calculated 
from the first two hundred prime numbers and the as- 
sociated semi-classical potential from equation ([7]). It is 
apparent that the inversion potential oscillates around 
the semi-classical potential except close to the edge of 
the potential. Similarly, one may solve the Schrodinger 
equation with the semi-classical potential and with the 
inversion potential obtained above, comparing how well 



4 



(a) ' ' ' f(b) ' ' ' ' 




Distance, x Distance, x 



FIG. 1: The reference potential, Vo{x) — x (dashed line), and the inversion potentials (solid lines), V{x), are shown using (a) 
one (b) two and (c) five energy eigenvalues indicated with horizontal dashed lines. The lower figure depicts the reference and 
inversion potential derived using the first fifty energy eigenvalues; the inset illustrates the difference between Vo{x) and V{x). 




Distance, x Distance, x 

FIG. 2: Figures show the reference potential, Vo{x) = x"^ , (dashed line) and the inversion potential, V{x), (solid line) using 
(a) one (b) two and (b) five energy eigenvalues which are indicated with horizontal dashed lines. In the lower panel, the same 
quantities are presented with the first one hundred energy eigenvalues utilised. The inset depicts Voix) — V{x). 



they reconstruct the original set of eigenvalues. Table 
im comprises a selection from the original set of eigenval- 
ues (eo,n)) labeled by n, and the energy eigenvalues of the 
semi-classical (esc.n) and inversion potential (e„) with the 
same quantum number. This also served as a numerical 
check of our method. The semi-classical energy eigenval- 
ues capture the trend, but - as expected - those derived 
from the inversion potentials are much better approxi- 
mations to the exact eigenvalues. One may notice that 
the agreement between eo,n and e„ is consistently much 
better for even values of n. The same tendency can also 
be seen in Tablc|T]for the triangular and harmonic oscilla- 
tor potential, although on an order of magnitude smaller 
scale. The underlying reason for this effect is the subject 
of ongoing investigation. 



C. Inversion potential for C,{s) zeros 



Similarly to the semi-classical approximation derived 
for the prime numbers, one may calculate a potential 
corresponding to the Ricmann C,[s) zeros using the fact 
that the number of zeros [l^ 



M{E)^—E\n{E)- ' E+- + 0{HE)) (9) 



Calculating the average density of the C(s) zeros from 
the expression above restricts the choice of the otherwise 
arbitrary reference energy level to i?o > Stt. Inserting the 
density into the Wentzel-Kramers-Brillouin semi-classical 




Distance, x 



FIG. 3: The main figure depicts the semi-classical potential 
(dashed line), Vo{x), and the inversion potential (solid line), 
V{x), derived from using Marchenko's method with the first 
two hundred prime numbers. The inset shows the difi'erence 
of Vo{x) and V{x). 




Distance, x 

FIG. 4: Main figure shows the semi-classical potential (dashed 
line), Vo(a;), and the inversion potential (solid line), V{x), 
derived from the inverse scattering method using the first two 
hundred energy eigenvalues. The inset depicts the difi'erence 
ofVo{x) and V{x). 



quantisation formula, we obtain (see 2.727.5 in [l8| ) 



x{V) = - 

TT 



VV - Eo In 



2^ 



\Vv~Vv-Eo 



(10) 



which is identical to that given by Wu and Sprung [l9|. 
The structure of the semi-classical potential close to the 
origin depends on the choice of the reference energy level. 
If Eq > 2-K then V(x) — Eq ^ , but in case of E^ — 
2TT the potential grows as V{x) — Eq ^ x^l'^ . As \x\ 
approaches infinity the potential becomes independent 
of Eq and expression (jTU]) can be inverted to obtain the 
following asymptotic 



V{x) 



w, 4^ 



(11) 



where W{z) denotes the Lambert-W function. 

Applying the formulae (jSal c) we calculated a number 
of potentials supporting the first n zeros of the C(s) as 
energy eigenvalues, utilising a tabulated form of the low- 
lying zeros [i^]. As an example, we have plotted the 
potential for N = 200 in Figure |4] and in Table [ll] one 
can compare how well the energy eigenvalues of the in- 
version potential coincide with the original eigenvalues, 
i.e. the zeros of the Riemann C(s) function. In this case 
the agreement is even better than it was for the prime 
numbers. This can be explained by the much slower in- 
crease of the potential than that for prime numbers, as 
X approaches infinity. Similar effects are seen for the two 
pedagogical examples in Table IH The energy eigenval- 
ues for the triangular potential are, at least, an order of 
magnitude more accurate than those for the harmonic 
potential. 



III. COMPARISON WITH EARLIER RESULTS 

Both variational and dressing-transformation tech- 
niques, have already been applied to construct quan- 
tum mechanical potentials for which the energy eigen- 
values are either the zeros of the Riemann C(s) function 



[13l . Il9l . l21j . or the prime numbers [2]J. The common 
feature of these methods is that the potential is built up 
in recursion by incorporating more and more eigenvalues 
into the spectrum. 

Previous works [l^ [2l| estimated the box-counting di- 
mension of the potentials belonging to the prime numbers 
to be 1.8, and for the Riemann C,(s) zeros to be 1.5, where 
the number of eigenvalues used ranged from 100 to 32000. 
Our measurements broadly support these values (see -Dq 
in Figure O. 

In order to reproduce these findings we treat the graph 
of the potential as a signal. To measuring the fractal 
dimension we de-trend the signal, i.e. subtract the ac- 
tual inversion potential from the semi-classical potential, 
f (x) — V{x) — V^^{x). Moreover, we limit ourselves to the 
spatial range of [0, 10] to eliminate any boundary effect 
arising from the fact that both Marchenko's method and 
the dressing-transformation produce a potential with a 
constant asymptotic value for large spatial coordinates. 

We have measured the Renyi-dimension of the poten- 
tials d, [12], defined as 



a — 1 



lim 

£^0 + 



In 



In (6) 



(12) 



where jti is the probability that the discrete random 
variable X falls into a box centered at Xi with side e. 
This probability can be approximated using the rela- 
tive frequencies obtained by dividing the embedding two- 



dimensional [x, V) space into a finite number of bins and 
counted how many times the potential takes its value 
in the given box. Contrary to the general box-counting 
method, which treats every box equally regardless of how 
many points of the fractal a given box contains, if a > 
in ()12p then boxes with higher relative frequencies will 
dominate the summand, therefore determining Da{X). 
On the other hand, if a < then the formula weighs the 
less dense boxes more and measures their scaling proper- 
ties. As a special case, a = associates equal weights to 
every box, and therefore Dq should reproduce the box- 
counting dimension. It can also be shown [9| that for 
a —^ 1 the numerator including the pre-factor, converges 
to the Shannon-entropy, defined to be — X^i-P* (P*)- 
Concluding, by calculating the generalised Renyi dimen- 
sion one can "scan" the fractal by its "density" , and mea- 
sure its heterogeneity. If Da varies with a then the fractal 
is a multi-fractal, since its subsets scale differently. 

One may derive two statistics for ^ based on p^ : (a) 
using a two-dimensional grid and approximating the two- 
dimensional conditional probability distribution with the 
relative frequencies of the individual boxes, or (b) calcu- 
lating the generalised dimension for the marginal proba- 
bility distribution of ^. 

In Figure [5] both set of statistics are demon- 
strated showing the generalised dimension for the two- 
dimensional probability distribution and the insets de- 
picting Da for the marginal distribution. The box- 
counting dimensions, Dq, are also indicated on the 
graphs. Although both curves in the main figures have 
an overall "S" shape, their structure is different. In the 
case of the potential generated from the prime numbers, 
both the positive and negative a branches of the curve 
show monotonic decrease towards the right. However, 
for the potential designed from the zeros of the Riemann 
^{s) function, the negative a branch of the curve remains 
nearly constant. 

These results suggest that the potentials calculated for 
the prime numbers and for the zeros of the Riemann 
C(s) function are indeed multi- fractals [1^. The steep 
middle par t of the curves also explains why earlier stud- 
ies [H, In, [U, [Hi differed in the box-dimension. The 
number of incorporated energy eigenvalues strongly in- 
fluence the conditional probability associated to one box 
and eventually shift Dq. 

Finally, we mention another conjectured property of 
the quantum system supposed to possess the zeros of 
the Riemann C(s) function as energy eigenvalues, namely 
that Ti. may violate the time-reversal symmetry. 

Similarly to the prime counting function, tt{x), one 
may define a function, Af(t), which counts the zeros 
of the Riemann (^(s) function, i.e. a function jumping 
by unit y w henever t passes over of the zeros, t„. It is 
proven [1^ that the function Af can be decomposed into 
a smooth and fluctuation part: Af{t) = N{t) +Afosc{t), 
where N{t) has been given in © explicitly. The fluctuat- 
ing term has remarkably similar structure to Gutzwiller's 
trace formula 25] giving the density of states of a quan- 
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a 

FIG. 5: The Renyi-dimension, Da, is shown for the potentials 
belonging to (a) the prime numbers and (b) the Riemann (^(s) 
zeros. Insets on the left show Da for the marginal distribution 
(see text) while those on the right depict the de-trended data. 
Here, in the dressing-transformation, we used the first 10^ 
eigenvalues and spatial step-size h = 10"^ for both the prime 
numbers and the zeros of ({s). 



tum system. The comparison of the two formulae led 
to the hypothesis (26j that a quantum system with the 
zeros of the Riemann (^(s) function as energy eigenval- 
ues, does not possess time-reversal symmetry. The ap- 
proach presented in this paper creates a symmetric, one- 
dimensional, although multi-fractal potential, V{x), for 
which the corresponding energy eigenvalues coincide with 
the first n non-trivial zeros of the Riemann (^(s) function. 
The reflection symmetry of the potential, V(x) = V{—x), 
guarantees time-reversal symmetry of the Hamiltonian, 
H = /2m + V{x). This result, therefore, allows us to 
assume the existence of a quantum system having the C 
zeros as energy eigenvalue and obeying the time-reversal 
symmetry simultaneously. 
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IV. CONCLUSION 

In the present paper we used Marchenko's method, 
one of the inverse scattering methods, to construct one- 
dimensional, symmetric quantum potentials, the energy 
eigenvalues of which coincide with either the prime num- 
bers or the zeros of the Riemann ({s) function. We have 
demonstrated the accuracy and usefulness of this method 
on two pedagogical examples, the triangular and har- 
monic potentials. For both cases we showed the recon- 
structed potentials and calculated the energy eigenval- 
ues, which agreed with the predescribed values very well. 
Later, we applied the same technique and numerically 
calculated potentials for the prime numbers and zeros 
of C{s). We found that the outcome of the Marchenko 
method is identical to that of the dressing-transform used 
previously. At the present stage, the latter method is nu- 
merically preferable to the Marchenko's method. Using 
the dressing-transform, we created potentials, to high ac- 
curacy, from the first one-hundred thousand prime num- 
bers and also for the same number of zeros of the Rie- 
mann C(s) function. Looking at the graphs of these po- 
tentials as signals one can analyse their statistical prop- 
erties. After de-trending these signals we calculated the 
Renyi-dimension, which is a generalised fractal dimen- 
sion. Our results suggest that inversion potentials are 
multi-fractals for both the prime numbers and for the 
zeros of C(s). The specific values of the generalised di- 
mension for the prime numbers Dq = 1.808, and for the 



C(s) zeros, Dq = 1.458, agree well with [2l[. 

Even though Marchenko's method is not yet able to 
compete with the dressing-transform in the number of 
eigenvalues incorporated into the potential, it gives ex- 
plicit formulae for how one can build up such potentials 
without recursion. Looking at formula (jSap one may see 
that the determinant of the matrix I + C is a polyno- 
mial of the entries, i.e. of exponential functions. Taking 
the natural logarithm and differentiating twice with re- 
spect to the spatial variable will not change the fact the 
potential is a rational function of a finite number of ex- 
ponentials. This fact is more apparent in formula ([6]). 
This gives us hope to be able to explore analytically the 
properties of a quantum system which possess the zeros 
of the Riemann C(s) function as energy eigenvalues. 

Simplification of formulae (|5al -c) is under investigation. 
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